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Abstract 

We compute the small time asymptotic of the fundamental solution of Hormander’s type 
hypoelliptic operators with drift, at a stationary point, xq, of the drift field. We show that 
the order of the asymptotic depends on the controllability of an associated control problem 
and of its approximating system. If the control problem of the approximating system is 
controllable at Xq, then so is also the original control problem, and in this case we show that 
the fundamental solution blows up as where iV is a number determined by the Lie 

algebra at a:o of the fields, that define the hypoelliptic operator. 


1 Introduction 

In the following paper we will study the small time asymptotic on the diagonal of the fundamental 
solution to an hypoelliptic operator with drift. Let M be a n dimensional smooth manifold and 
let /o, fi,...,fk be smooth vector fields on M, where /o is usually called drift field. We will 
consider the following operator on X M) 

k 

e C-(M+X M). (1) 

2=1 

We will assume that the vector fields fo, fi, ■ ■ ■, fk satisfy Hormander’s condition of hypoel- 
lipticity (see m), that is 

Lie,{^=K xT,M VxGM. (2) 

Under this hypothesis the operator in ([1]) is hypoelliptic and, as shown for example in [9] and 
[8], it admits a well defined smooth fundamental solution, p{t,x,y), for small time, that is given 
by the probability density of the stochastic process at time t, starting at x at time 0 and 
driven by the stochastic differential equation 

k 

df,t = k{ft)dt + X fiiit) o dwi 
2=1 

where w = {wi,... ,Wk) denotes a /c-dimensional Brownian motion and o denotes the integration 
in the Stratonovich sense. 
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Given a stationary point, xq, of the drift field /q, i.e. fo{xo) = 0, the goal of this paper is to 
study the behaviour of the heat kernel, p{t,x,y), for small time, in xq. 0 

This problem has already been studied by many authors in the past and recent period. In 
particular, if the vector fields /i,..., /^ are Lie bracket generating, i.e. they satisfy condition ([2]) 
without the contribution of /o, many results are already known. Indeed under this hypothesis 
it is well defined a distance function d{x,y) determined by the fields fi,..., fk- Then it was 
proved by Leandre in |10| and m that 

lim t logp{t,x,y) = -^(f{x,y) 
t^O 2 


which generalizes a result by Varadhan m in the elliptic case. If moreover /q = 0, then if 
X ^ y and x and y are not conjugate along any minimal geodesic, then for some constant 
C = C{x,y) e M 


P{t,x,y) 


C + Oit) 

-2t 

fn/2 ■ 


as proved by Ben Arous |3] and Barilari, Boscain and Neel |3]. Here we see in particular that 
the order of the asymptotic in t is determined by the dimension of the manifold. On the other 
hand if x = y and fo = 0, Ben Arous and Leandre showed in [5] and that 


p{t,x,x) 


C + 0{Vi) 

tQ/2 


(3) 


where Q is the homogeneous dimension of the manifold and C is a constant depending on x. 
These two examples show that the order of degeneracy of the fundamental solution for small 
time can reveal interesting geometric properties and depends on the structure of the manifold. 

In this paper we are going to determine the order of the asymptotic for small t and we will 
see that if /o is not identically zero the asymptotic Q on the diagonal at a stationary point, 
Xq, of the drift field could be more general, as it was already pointed out by Ben Arous and 
Leandre in |6]. The proof relies on a homogeneity property of the operator in ([1]) under the 
dilations. Indeed to understand better the behaviour of the fundamental solution p(t, xo,xo), 
we will dilate the space around xq and rescale the time. We will see that, under a suitable 
dilation, this procedure splits the operator ^ ~ fo ~ ^ Yli=i ff ^ principal operator 



2=1 


(4) 


perturbed by a small operator that goes to zero as the space becomes larger. Here the fields 
/o, /i,..., /fc are nilpotent approximations of the original fields fo, fi, ■ ■ ■, fk around the point 

^The assumption of fo{xo) = 0 is important, but can be avoided in some cases. For example, one could be 
interested in studying the differential operator 

^ - /o - ^ 

i=l 

that can be treated as an operator of the form (B 

i = l 

once we take fo~fo + ^ X]iLi(hiv^/i)/i. The results of this paper apply also to an operator of this form and we 
will show the small time asymptotic of the fundamental solution on the diagonal at a point xo, where fo{xo) — 0 
(not necessarily fo{xo) = 0). Details for this kind of operator are written in Remark 15.51 
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xq. In this procedure it is very important to find a good coordinate system around xq, so that 
the dilation of the space produces an hypoelliptic principal operator (j?]) . This is done by noting 
that the field /o acts in the operator just one time, while the fields fi, ■ ■ ■, fk are applied twice. 
So we will give to the drift field a double weight with respect to the others. We explain this 
procedure in detail in Section [2j 

We will prove that the asymptotic of the fundamental solution depends on the controllability 
of some associated control problems: consider the control problem 

k 

x{t) = fo{xit)) + Mt)Mx{t)) (5) 

i=l 

where u = {ui,... ,Uk) G L°°(]R^;M^) are the controls. If the control problem ([5]) is not 
controllable around the point xq, i.e. from xq we can not reach a neighborhood of xq using 
curves described by the control problem, then 

p{t, Xo, Xq) =0 Vt > 0. 


We will consider also the control problem induced by the approximating system 


k 

x{t) = foix{t)) + ^ Ui{t)fi{x{t)). 
i=l 


( 6 ) 


If it is controllable, then so is also the original one ([5]), and, if we denote by qQ{t,x,y) the 
fundamental solution of the approximating hypoelliptic operator in ([1]), we prove in Theorem 
15.41 that 

/, X qo{'i-,xo,xo) +0{t) 

p{t, Xo, xo) = --, (7) 


where the factor comes from the change of coordinates given by the dilations and it is 
determined by the order of homogeneity under the dilations of a volume form like dxi A... A dxn 
around xq. 

The number N is an integer depending on the structure of the Lie algebra generated by the 
helds /o, /i,..., /fc in Xq. It is computed in detail in Section 01 but to anticipate here how it is 
found, we give a brief explanation of it: for j G N, let 


dj := dim span,^^, {[/ii, [• • •, [/ij.i, /ij • • ■]] such that #{4 > 0} + 2#{4 = 0} < j} 

be the dimension of the space generated by all the Lie brackets at xq of at most j vector fields 
/o, /i,..., /fc, where the field /q is counted two times. Then the number N is defined as 

i>i 

For an operator coming from a sub-Rimemannian manifold with no drift field, N is equal to the 
homegeneous dimension Q of the manifold. If the drift field is not zero, then N could be more 
general. 

In the intermediate case in which the approximating control problem ([6]) is not controllable in 
Xq, but the original control problem Q is still controllable, then the behaviour of the asymptotic 
can be more general. It could vanish, or blow up faster than and even exponentially fast, 

as it was already pointed out by Ben Arous and Leandre in [6]. 

The paper is organized as follows. We will begin by describing in details in Section [2] the 
homogeneity properties of the operator in ([1]). In particular we will derive the conditions that 
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the dilations have to satisfy in order to produce the right split of the dilated operator, into a 
hypoelliptic principal part, perturbed by a small operator. We will also give a brief introduction 
into Duhamel’s formula in Subsection 12.11 which will be an important tool to study the perturbed 
operator. In Section [3] we will introduce the coordinates that give the right dilations of the space. 
This coordinates are dehned from a hltration of the tangent space to xq determined by the helds 
/o, /i,..., /fc and give rise to a graded structure around xq, which defines an anisotropic dilation. 
In Section [I] we will define the nilpotent approximation, that determines the principal part of 
the vector fields /o,/i,...,/fc, which define the principal operator (j?]). We will compute also 
the integer N appearing in the asymptotic dZD, that derives from the change of the volume 
form under the dilations. In Section [5] we will prove the asymptotic ([7]), by using the tools 
introduced in the previous sections. In the following Section [6] we will focus our study on the 
operator derived from the nilpotent approximating system and its associated control problem. 
By proving a modihcation of Stroock and Varadhan’s support theorem, we will give a necessary 
and sufficient condition for the positivity of the fundamental solution of the principal operator, 
that is based on the controllability of the approximating control system. Finally we will end the 
paper with Section [71 where we show a series of examples, illustrating how this formula recovers 
in particular the known results recalled in the introduction. 

2 The fundamental solution and its behavior under the action of 
a dilation 

Let M be an orientable n-dimensional manifold (without boundary) and let // be a volume 
form on M. Given fo, fi, ■ ■ ■, fk smooth vector helds on M we consider the following partial 
differential operator: 


V<^eG-(MxM). (8) 

i=l 

We will call /o the drift held and we will denote by L the operator /o + ^ ff- Let us 

recall the dehnition of fundamental solution. 

Definition 2.1. The fundamental solution of an operator — C over M x M with respect to 
the volume p, is a function p{t,x,y) G X M X M) such that 

• for every fixed y G M, it holds x, y) = Cxp{t, x, y), where the operator C, acts on the 

X variable; 

• for any ipQ G C^{M), we have 

lim / p{t,x,y)(fo{y)Tiy) = To{x) uniformly in t. 

* xO Jm 

In other words, if we want to solve the partial differential equation ^ = C(p with initial condition 
99 ( 0 , x) = (Pq{x), the fundamental solution allows to reconstruct ip by convolution of po with 
p{t,x,y). 

Remark 2.2. The fundamental solution depends on the given volume p in the following way: 
let p and v be two volume forms on M, and let ghe & smooth function such that v = e^p. Let 
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Pfj_ and pi, denote the fundamental solutions of ^ £ with respect to p and respectively. Then 

for every initial condition ipQ S C^{M), the solution ip{t,x) of 


I £(0,x) = ipo{x) 


is given by 


ip{t,x)= Pf,{t,x,y)ipo{y)p{y) 

J M 

= Pu{t,x,y)ipo{y)i^{y)= pu{t,x,y)ipo{y)e^^y'> p{y), 

J M J M 


where the equalities follow since the solution is unique for smooth vector fields, 
arbitrary, we have 


p^{t, X, y) = X, y) Vt > 0, Vx, y G M. 


Since ipQ is 


From the point of view of the asymptotic of the fundamental solution on the diagonal, it 
follows that the two asymptotics are the same for both volume forms up to a multiplicative 
constant 7 ^ 0 depending on the relation between the two volume forms and on the point 

where we compute the asymptotic. For the study of the small time asymptotic on the diagonal, 
we will then suppose p = dxi A ... A dxn near the point xq. 

The existence of a fundamental solution for an operator ^ £ is not always guaranteed, 

but in 1967 Hormander published an important paper on the properties of an operator like the 
one in ([8|). An operator £ is called hypoelliptic if whenever, for a function tp, Cp G C°°(U) on 
an open set U, then cp G C°°(U). In [ 8 ] Hormander proved a condition of hypoellipticity for 
operators of the form ([ 8 |), namely if 

Lie,{^-/o,/i,...,/fc} =MxT,M VxGM, 

then the operator ^ — L is hypoelliptic. Throughout this paper we will always assume that the 
fields /o) /i) ■ ■ ■, /fc satisfy Hormander condition. 

Under this assumption of hypoellipticity, the operator in ([ 8 ]) admits a fundamental solution, 
p{t,x,y) (see the backward Kolmogorov equation and [9], Chapter 4.6 for the study on R” and 
Chapter 5 for the equation on a manifold). It is given by the probability density function of the 
process that satisfies the following stochastic differential equation in Stratonovich form 

I d^t = fo{Ct)dt + Yli=i M^t) o dwi, 

\ (o = x 


where is a process on M and {wi,... ,Wk) denotes a ^-dimensional Brownian motion. More 
explicitly, the fundamental solution of dH) with respect to a volume form p on M is the function 
p{t, X, y) such that for every measurable set A C M, the probability of the process to be in A 
at time t is given by 

^[6 e = 3 :] = / P{t,x,y)p{y). 

Ja 

Moreover, for every initial condition ipQ{x)^ the corresponding solution of the differential equation 
is 

p{t,x)= / p{t,x,y)po{y)p-iy) = = x], 

J M 
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that is the expectation value of (po{^t) at time t, knowing that = x. 

Let xo be a fixed point of the drift held, /o(xo) = 0, our goal in this paper is to understand 
the small time behaviour of the fundamental solution at the point xq- Our method will proceed 
as follows: we will dilate the space around the point xq, and we will rescale the time accordingly, 
so that the space around xq will be magnihed and we can study better the behaviour of the 
solution at xq- 

Definition 2.3. Let (U,x) be any coordinate neighborhood of xq and let {wi,... ,Wn) be positive 
integers, that we will call weights of the coordinates (xi,... ,Xn)- 

For €> 0 we define the dilation, 6 ^, of order e and weights [wi,... ,Wn) around xq as the 
function 5^ : 17 —)• [/, such that 


(5(xi,...,Xn) := Vx = (xi, ..., x„) e C/. (11) 

For e> 1 we understand the dilation defined actually on a smaller neighborhood of xq in U. 


Under the action of the dilations 5^, the coordinate functions and the coordinate vector helds 
behave as 

-| O 

XiO — = — 7 ^— VI < z < n. (12) 

-^*dxi e'^idxi - - ^ ^ 


Let the volume p, be represented in the coordinate neighborhood (U, x) by /r = dxi A ... A dxn- 
By Remark 12.21 this assumption is not restrictive for the study of the asymptotic along the 
diagonal. Then the volume form /i changes under the action of the dilation as 


5* (dxi A ... A dXn) = e^»=i '^^dxi A ... A dXn- 


(13) 


When we apply a dilation to the space around xq and we rescale the time variable, also the 
fundamental solution is changed accordingly, as it is proved in the following proposition. 

Proposition 2.4. Let [U, x) be a coordinate neighborhood around the point xq and let pL be a 
volume form on M such that pL = dxi A ... A dxn in U. For weights (rci,..., Wn) and 0 < e < 1 
consider the dilation 5^ ■ U —)• U. Let p{t,x,y) be the fundamental solution of the operator in 
([8|) with respect to the volume pi. Then the fundamental solution on U of the operator 

|-cU/../o+i|;(i./../.)7 (M) 

is the function 

qe{t, X, y) := e^^=i p(e“t, 5e{x), 6 ^{y)) Vx, y G U. (15) 

Here denotes the pushforward of a vector field X under the action of 6 ^ and a is a real 

positive number. 

Remark 2.5. The coefficient of normalization eEi=i"'7 that we have used to dehne q^, is 
necessary in order to reconstruct all the solutions of the differential operator, by convolution 
with q^. This coefhcient appears as soon as we make a change of coordinates in the integral 
of the convolution. Moreover, we will see that from this coefhcient we hnd the order of the 
asymptotic of the fundamental solution for small time. 

Proof. First of all notice that the dilation 61 /,, : d^{U) —)■ U can be dehned only on the smaller 
neighborhood S^{U) of U. Then the helds Si/^^:fi are vector helds just on the coordinate neigh¬ 
borhood U. 
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Next let us prove the first property of the fundamental solution, i.e. that the function is 
a solution of the operator in (I14p . For convenience, we call ip the dilation from M"*" x Si/^U to 
M"*" X U defined by 'ip{t,x) := (e“t, (fe(x)). Then the function can be written as q^{t,x,y) = 
p{ip{t,x),5e{y)) and the operator in (fTTll is 


V’* “ V’* Vo 




where we are using a little abuse of notation, by considering /q, ..., /fc as vector fields defined 
on the product space M"*" x U. Recall the definition of pushforward of a vector field X under 
the action of a diffeomorphism (p: for every function g we have Px*{X)if)\ip(^x) ~ ^xif ° ^\x)- 
Then we compute for fixed y ^ U 



qeit,x,y) 




{p{'ip{t,x),5^y)) 


P\ 


dt 


P\{ip(t,x),Sey) 


{■tj)(t,x),&ey) 


{ilj(t,x),5^y) 
e^'-l M(qp{t,x),y) P\{ilj(t,x),5^y) 


where the last equality follows since p is the fundamental solution of the operator ^ — L. 
Applying the same computations to the vector fields for i = 0,...,k, we find that q^ 

satisfies 

d 1 ^ 

'ipP^-^Qe{t,x,y) = ipp^fo{qe(t,x,y)) + 2'^ V^Vi) {qe{t,x,y)) 

i=l 

and hence q^ is a solution for the operator in (1141) . 

Let us prove the second property of a fundamental solution. Here we will see that the 
constant of normalization in the definition of q^ is exactly the parameter that we need 

in order to construct the other solutions of the partial differential equation by convolution with 
the fundamental solution. Indeed let us prove that for any ipQ £ C^(U), it holds 


lim / q^{t,x,y)po{y)p{y) = Po{x). 
i\0Ju 

This follows by a change of variable and the same property valid for the fundamental solution 
p{t,x,y): 

lim / q^{t,x,y)ipo{y)p{y) = lim / e^'^=^'^^p{eH,5eX,5ey)po{y)p{y) 

= lim f e^i=i^y{e'^t,deX,6ey)po{y)p{y) 

We can integrate on M, because, since (fQ has compact support in U, we can extend (po to be 
zero outside U. Now let us do a change of variable with 2 : = S^y- As computed in (1131) . the 
volume is transformed as p{z) = e^*=iThen 

1 ™ / q,{t,x,y)po{y)p{y) = liin / (f.x, z)(y?o o Si/t){z) 

t\oJu t\oJM 

= {^ 0 ° ^l/e){^eX) = Po{x) 

where the second equality follows because p is a fundamental solution. □ 
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Let us investigate better the behaviour of the fields /j under the action of the dilations. Let 

(i) 

us write every component of the fields /j in a Taylor expansion centered in xq = 0 for x in 
the coordinate neighborhood U 


f 

I 


U) 




^ 8 #( 0 ) 

l=l 


dx 


d 


By the properties of a dilation acting on the coordinate functions and on the coordinate vector 
fields, (1121) . when we apply a dilation to the vector fields fi, every component has a different 
degree with respect to e. Depending on the value of the weights (rci,..., and on the 
coefficients of the Taylor expansion of the fields fi, for every i = 0,... ,k, there exist an integer 
Oj and a principal vector field fi such that 


hji = — fi + o 


1 


where fi contains the components of every fi^^-£^ that is homogeneous of degree —a* with 
respect to the dilations. Applying this formula to the dilated operator in ([HI) , we find that the 
operator L rescales as 


^ {^l/e*fi) 


2 = 1 


1 - 1 ^ 1 

+ 2 X] ° 
2=1 



where a := max{aO) 2 q;i, ..., 2ak}. 

The main task in our study will be to find suitable coordinates and good weights Wi, so that 
all the principal parts of the vector field /o and of /f,..., rescale with the same degree under 
the dilations. In this way the principal part in the dilated operator is homogeneous of order —a 
and if we choose a = a, the operator L in (fHl) can be written as 


/ , k \ k ^ 

( ^l/e*fo + X ^ {^l/e*fi) 1 = fo + ^ fi + (16) 

V ^ i=l / ^ i=l 

where is a differential operator that goes to 0 as e goes to zero. As proved in Proposition 12.41 
its fundamental solution on U is given by 

qe{t, X, y) := S^x, 5^y). (17) 


To an operator like the one in (I16D we can apply Duhamel’s formula, that gives the asymptotic 
of the fundamental solution as a perturbation of the asymptotic of the fundamental solution of 
the principal operator. 


2.1 Duhamel’s formula 


In this section we recall briefly a famous formula, called Duhamel’s formula, which allows to find 
the asymptotic of the fundamental solution of a perturbed operator, once we have the explicit 
fundamental solution of the principal part of the operator. This method is presented, among 
others, in Chapter 3 of |12| and in [2]. 

Let C be an operator on a Hilbert space with fundamental solution p{t, x, y) and let us define 
the following operator on the Hilbert space 


e 


tc 



p{t,x,y)ip{y)dy. 



By the properties of the fundamental solution this is an heat operator i.e. an operator such 
that 


—-— = Ce^'-'ip and lim 

dt ^ t^o ^ 


Suppose that C can be decomposed in a sum, 

C, = ds) + X , 

of a principal part, /Iq) a-^d a perturbation, X, and assume that Cq has a well dehned heat 
operator Then Duhamel’s formula allows to reconstruct the heat operator of £ by a 

perturbation of the heat operator of £o (s 66 Chapter 3 of m for a proof), namely 


^tC _ ^tCo _|_ f g(t s)Cp^^sCo^g _ ^tCo _|_ ^tC ^ 

Jo 


(18) 


where with * we denote the convolution operator between two operators, A{t) and B{t), on the 
Hilbert space: 


{A*B){t)= f A{t — s)B{s)ds 

Jo 


Let a{t,x,y) and h{t,x,y) be the heat kernels of A{t) and B{t) respectively and let X be an 
operator. Then the kernel of {A * XB){t) is obtained as follows: let 99 be a function in the 
Hilbert space, then 


[{A*XB)it)y,] {x) = 


[ A{t — s)XB{s)ds ip 

Jo 

f A{t-s) X f b{s,-,y)(p{y)dy 
Jo L Jm 


ds 


n a{t - s,x,z)xA / b{s,z,y)ip{y)dy] dzds 
4 \J M / 

i a{t - s,x,z)Xzb{s,z,y))dzds] (p{y)dy 

Jm \Jo Jm J 


so the heat kernel of {A* XB){t) is 

{a * Xb){t,x,y) = / / a{s,x, z)Xzb{t — s, z,y)dzds. 

Jo Jm 

From (fT 8 |) we can now derive an approximation of the heat kernel p{t, x, y) of the perturbed 
operator C, by means of the heat kernel pQ{t,x,y) of the principal operator: 


p{t, X, y) = po{t, X, y) + {p* Xpo){t, x, y). 


(19) 


We will make use of this formula in the next section to find the asymptotic expansion of the 
fundamental solution of ©• 


2.2 The perturbative method 

We can apply Duhamel’s formula to the operator in (11611 . Indeed it is the sum of a principal 
operator 

lq = fo+ 

i=l 
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perturbed by a small operator If we find good coordinates and weights, so that ^ — Lq 
admits a fundamental solution qQ{t,x,y), then by Duhamel’s formula (1191) the asymptotic of the 
fundamental solution is 


qe{t, X, y) = qo{t, x, y) + {q^ * X^qo){t, x, y). 


( 20 ) 


From here we will immediately conclude our study, indeed let us choose in the definition of 
dlZD, X = Xo, let e go to zero as e = and fix the time variable i = 1, then we get 


p(t,xo,xo) 


ge(l,xo,xo) 


1 




(go(l,3:o,3:o) + * <^ego)(l, a;o)) • 


In conclusion, in the choice of the coordinates {U,x) and the weights {wi,... ,Wn) it will be 
important that the principal parts, fo, fi, ■ ■ ■, fk, of the vector fields make homogeneous the 
principal part of the dilated operator and, moreover, that they satisfy Hormander condition, so 
that it is guaranteed the existence of a fundamental solution go of the principal operator Lq. 
Finally we will need to check that the remainder term (g^ * df(;go)(l, xq,xq) in the asymptotic of 
p goes to zero, as t (and therefore e) goes to zero. 


3 Graded structure induced by a filtration 

In this section we will introduce some notation and recall the definition of local graded structure 
of a manifold, induced by a filtration. This terminology will be essential to find the right 
coordinate to rescale the differential operator L and to compute the order of the asymptotic of 
the fundamental solution. We will constantly refer to Bianchini and Stefani’s paper [7]. 

3.1 Chart adapted to a filtration 

Let M be a n dimensional smooth manifold and /o, /i, • • •, /fe be smooth vector fields on M, that 
satisfy Hormander condition ([2]). We will consider the hypoelliptic operator on M x M defined 
as 

h-h-hi 

2 = 1 

As noticed in the previous section the role played by the drift field /o and the other vector 
fields, {/i,..., fk}, in the sum of squares, is different, and in particular the fields {/i,... , fk} 
are applied twice as many times as the drift field is. For this reason we want to treat differently 
the two kinds of fields, by giving to them two different weights. 

Let LieAl be the Lie algebra on M generated by a set {Xo,Xi,... ,Xk} of noncommutative 
indeterminates. 

Definition 3.1. For every bracket A in LieAl we denote by |A|j the number of times that the 
indeterminate Xi appears in the definition of A. We will call this number the length of A with 
respect to X^. 

For example, the bracket A = [Xq, [X 2 ,Xo]] has lengths |A|o = 2, |A|i = 0 and |A |2 = 1, and 
it has zero length with respect to any other indeterminate. 

By fixing a weight, fj, to every indeterminate AIq, ..., Xk be can define the weight of a bracket 

A. 
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Definition 3.2. Given a set of integers ... ,1^), we define the weight of a bracket A E 

LieX as 

k 

||A|| ■.= Y,W^ 

i=0 

and we set ||0|| = 0. 

In order to give different importance to the drift field, with respect to the other vector fields, 
in the following we will fix the integers to be 

Iq = 2 and li = ... = 1^ = 1. (21) 

This means that the indeterminate Xq will have weight 2, while the other indeterminates will 
have weight 1. For more complex Lie brackets, we have for example that the weight of the 
bracket considered before A = [Xq, [X 2 ,Xo]] is 

||A|| = 2-|A|o + 1-|A|2 = 5. 

By means of the weight of the indeterminates we introduce now a filtration of the Lie algebra 
spanned by /q, /i,..., /^ in the following way. For every bracket A in LieX we denote hy Aj the 
vector field on M obtained by replacing every indeterminate Xi with the corresponding field /j 
for 0 < i < k. Then we define an increasing filtration C = {Lj}j>o of Vec(M) by 

Li = span{Aj ; A G LieX, ||A|| < i}. (22) 

In other words Lj is the subalgebra of Vec(M) that contains all the vector fields obtained from 
a bracket of weight less then or equal to i. In particular, following our choice of weights, the 
first subspaces of the filtration are 


Lo = { 0 } 

Li = span{/i,... ,/fc} 

L 2 = span{fi,[fi, fj], fo : i, j = l,...,k} 

La = span{fi,[fi,fj],fo,[fi,[fj,fh]],[foJi]-i,j,h = l,...,k} 


Notice moreover that for every i,j > 0, the following properties hold 

• Lj C Lj_|_i 

• Ui>o = Lie{/o, /i, • • •, fk} and (J Lfix) = T^M, for every x G M, since by assump¬ 
tion the family {/o, ■ ■ ■ , /fc} satisfies the weak Hormander condition ([2]). 

When we evaluate C at the stationary point xq we get a stratification of the tangent space T^qM 
at Xo and this stratification will induce a very peculiar choice of coordinates around xq- 
Let 

di := dimLj(xo) Vi > 0. 

In particular, do = 0 and di < k. Moreover, by Hormander condition ([2]), there exists a smallest 
integer m such that L^ixo) = Tx^M. We will call this number the step of the filtration C at xq- 
The filtration C induces a particular choice of coordinates centered at xq, as proved by the 
following proposition. 
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Proposition 3.3 (Bianchiiii, Stefani [?])• There exists a chart {U,x) centered at xq such that 
for every 1 < j < m 

(i) Lj{xo) = span{^,..., 

(ii) Dxh{xo) = 0 for every differential operator D £ := {Zi---Zi with Zg £ Li^ and 

ii + ■ ■ ■ + il < j} and for every h > dj. 

Definition 3.4. We will call a chart that satisfies the properties of Promsition [Ql an adapted 
chart to the filtration C at xq. 

Since this kind of coordinates will reveal to be very important in our study, we give here the 
proof of the proposition, which relies upon the following Lemma: 

Lemma 3.5. Let m he the step of the filtration C, at xq and let j < m be an integer. Let 
(f £ C°°{M) be such that dxQ^p 7 ^ 0 and Z ■ (p{xq) = 0 for all Z £ Lj. Then there exists an open 
neighborhood U of xq and a function (p £ C°°{U) such that 

dxoP = dxo(p and D ■ (p{xo) = 0 VD £ = {Zi ■ ■ ■ Zi : Zg £ Li^,ii-\ - \-ii<j}. 

Proof. Let li,..., be vector fields on M such that they form a basis of Tx^M at xq and such 
that 


• {Li,..., Ydf\ are in Lj and form a basis of Lj at xq, for every i < j, 

• Yi ■ p{xo) = 0 for every i < n — 1, 

• Yn ■ p{xo) = 1 . 

We choose the chart y = (j/i,..., yn) as the local inverse of the map 

{yi, ...,yn)<-^ • • • e^^^^xo. 

Then the function <^ := is such that dx^p = dx^p- 

Let D = Zi ■■■ Zi £ A^ with Zg £ Li^ and L + • • • + z; < j and let us prove the second 
property required for the function p, by induction on 1. Since dx^p = dx^p and Z ■ p{xq) = 0 
for all Z £ Lj by hypothesis, the property is satisfied for I = \. For / > 1, since Zi £ Li^ we can 

write Zi(xo) = ^iYi(xo) for some Oj so 

dii 

D ■ p{xo) = ^ Oj (Zi-i ■ Yi ■ Zi _2 ■■■ Zi + [Yi, Zi_i] ■ Zi _2 ■ ■ ■ Zi) ■ p{xo) 

i=l 

The second component on the left side vanishes because, by the definition of the filtration, 
[Lj,L/i] £ Lij,.^, so we can apply on this component the induction hypothesis. By applying 
again the same commutation we have 

dii 

D ■ p{xo) = ^ aiZi_i ■■■Zi-Yi- p{xo). 
i=l 

Iterating the same procedure also to Zi_i,... ,Zi we can write D ■ p(xo) as a linear combination 
of elements of the type 

Yg-'-Yi^ -pixo), 

with 1 Y ifi Y dj < n for every h = 1,... ,1. Therefore we get D ■ p{xq) = D ■ ynixo) =0. □ 
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Proof of Proposition VJ.^A Let {U,x) be any chart centered at xq. We can get a chart with 
property (i) of the proposition by a linear change of coordinates. Let us still denote it by (U, x). 
For every i < n let j be such that dj < i < dj+i- Then the coordinate function Xi satisfies the 
hypothesis of the Lemma with respect to the integer j. By applying the Lemma to each function 
of the chart we get the statement. □ 

We will see an example of an adapted chart in Section 14.11 

3.2 Graded structure 

For 1 < f < m we define the integers 


ki 


di di—1) 


which indicate the number of new coordinates achieved with every new layer of the filtration at 
Xq. 

Definition 3.6. Let us denote a point x = (xi,... ,Xn) € U by the m-tuple ... ,x™) G 

0 0 ... 0 where each component x* := {xdi_^+i, ■ ■ ■, is a vector of length ki. 

Then for every e > 0 we define the anisotropic dilations around xq of factor e as 

Se{x) = ..., x™') := (ex^, e^x^,..., e^x”^). 


Notice that if e > 1, the map 5^ is actually defined only in a smaller neighborhood of xq, but we 
will still write the map 6^ U ^ U, understanding "locally defined". 


The dilations 5^ act on every coordinate function and on the coordinate vector fields with a 
different weight, namely 


Xj o 


= e*x 


j! 




^di—\ ^ j ^ di. 


For every 1 < f < n let Wj be the order of expansion of the coordinate Xj, that is wj := i if 
di-i < j < di. We call Wj the weight of the coordinate Xj. Then the dilation 5^ is a particular 
choice of the dilations defined in Definition 12.31 with respect to the coordinates induced by the 
filtration C and the weights {wi,..., Wn). 

Accordingly we define the weight of a monomial to be 


>V(x“ 


«i 


■xp") : = 


E 

i=i 


ajWj 


and the weight of a polynomial to be the greatest order of its monomials. Moreover we define 
the graded order, 0{g), of a function g G C°°{U) to be the smallest weight of the monomials 
that appear in any Taylor approximation of g at xq. 

For example, let n = 2 and suppose that xi has weight 1 and X 2 has weight 2. Then the 
polynomial X 1 X 2 — has weight 6, because the two monomials composing it are X 1 X 2 of 

weight 3 and the rest of weight 6. On the other hand, sin(xiX 2 ) = X 1 X 2 — —ho ((xi)^(x 2 )^) 

has graded order 3. 

We extend these definitions to differential operators. We say that a polynomial vector field 
Z is homogeneous of weight i if 


yV{Zip) = VV((/j) — i V monomial ip of weight >V(v9). 
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In other words Z subtracts weight i to every function. Then the weight of a polynomial vector 
field is the smallest weight of its homogeneous components. We define the graded order, 0{D), 
of a differential operator D by saying that 

0{D) < j if and only if 0{D ■ ip) > 0{ip) — j V polynomial p, 

that is D subtracts at most weight j from the functions. For example the graded order of a 
vector field like is 

O =Wh- ajWj 

Coming back to the previous example, the graded order of a field like sin(xiX 2 )g§^ is obtained 
as 1 — C)(sin(xiX 2 )) = —2. 

By means of the graded order we can give a generalization of the concept of Taylor approxi¬ 
mation of a function up to weight h. Namely, for any p £ C°°{U) and every integer h> 0, there 
is a unique polynomial of weight h such that 0{p — P(^h)) ^ 

Definition 3.7. The polynomial p^^p^ is called the graded approximation of weight h oi p and 
it is the sum of the polynomials of weight less then or equal to h in the formal Taylor expansion 
of p at xq. 

For every vector field V G Vec(C/) and each integer h <m there is a polynomial vector field 
of weight h such that 0{V — < h — 1. 

Definition 3.8. is called the graded approximation of weight h olV and it is the sum of 
the homogeneous vector fields of weight greater than or equal to h in the formal Taylor expansion 
of V at xq. 

Notice that, since V(/j) is a polynomial vector field, we can consider it as defined on the whole 
Euclidean space M"'. 

We will see in the next sections how to apply this graded structure in order to underline the 
most important properties of the operator in ([T|), concerning the small time asymptotic of its 
fundamental solution. 

4 Nilpotent approximation and the order of the dilations 

In this section we apply the graded structure, that we have just developed, to define a special 
class of vector fields, which approximate the original one /o,/i,...,/n and we will give an 
example to clarify the setting. Finally we will compute how the dilations change the volume 
form and we will introduce the order that will appear in the asymptotic of the heat kernel. 

4.1 Nilpotent approximation 

Let (x, U, w) be the graded structure around xq introduced in Section [3] and let fo, fi, ■ ■ ■, fk be 
the vector fields used to define the filtration of TxqM. Then as proved in [7] Theorem 3.1, for 
every / G Li, we have a bound on the graded order, namely 0{f) < i, where O is the graded 
order associated to the graded structure (x,U,w) defined in Section [3j 

Recall the integers Iq,Ii, ... ,lk introduced in Section [3] to define the filtration and denote by 
fi the graded approximation of weight /j of /j. In other words, fi has weight li and 0{fi — fi) < li. 
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Definition 4.1. The fields fo, fi, ■ ■ ■, fk called the nilpotent approximation of Jq, fi,..., fk- 

The fields fi, for 0 < i < k, are polynomials, so they can be defined on M”. 

We can describe more precisely the structure of the approximating fields ff. 

• /o contains the terms of weight 2; therefore every component /q of /o depends only linearly 
on the coordinates of weight Wj — 2 and more then linearly on the coordinates of less weight, 
but does not depend on the coordinates of weight greater then or equal to Wj — 1, that are 
Xh with h > dj- 2 - Moreover, since /o vanishes in xq, there are no constant terms. 

• fi contains the terms of weight 1 for i = therefore every component f- of /j 

depends only linearly on the coordinates of weight Wj — 1 and more then linearly on the 
coordinates of less weight, but do not depend on the coordinates of weight greater then or 
equal to Wj, that are Xh with h > dj-i. 


To make the construction more clear we end this subsection with an example, in which we 
will present the filtration in xq, the induced adapted chart and the graded structure, and we 
will find the related nilpotent approximation. 


Example 1. Let M = and let the number of controlled vector fields he k = 1. Define the 
vector fields 


d d 

■fi •= a- 

OX\ UX2 


and /o := sin(x^) 


A 

dx7 


and recall the choice of weights Iq = 2 and li = 1. The non vanishing Lie brackets that contribute 
to span the tangent space in any point are given by 


[/i,/o] = 2xicos(xi' 


A 

dxo 


[/iJ/i,/o]] = (2cos(x?) -4x?sin(x?)) 


A 

dxo 


Then Hormander assumption ([2|) holds in any point and the filtration defined in (|22p is equal to 


• Li= span{/i} 


• L 2 = span{/i,/o} 


• Ls = span{/i, /o, [/i, /o]} 

• L4 = span{/i, /o, [/i, /o], [/i, [/i, /o]]} 

In particular the filtration in xq is given by 

d 

Li(xo) = L 2 {xo) = Ls{xo) = span{-—} and L^ixo) = (23) 

UXl 

and the dimensions are: di = d 2 = d^ = 1, = 2. 

Let us find an adapted chart to the filtration at xq. As one can easily see, the coordinates 
(xi,X2) are not adapted, since (/i)^ • X 2 {xo) = 1 7^ 0 and the second property of the adapted 
chart then fails. Following the constructive proof of Lemma 13.51 one can find that the new 
coordinates (^1,^2) defined by 

3,2 

yi = xi - ^ + X2 

3,2 

2/2 = + X2 

give an adapted chart at xq. In this coordinates the two vector fields are written as 
/. = ^ and /„ = (A + 
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For e > 0 the dilations defined in Definition 13.51 are 


: (2/i, 2/2) (eyi,e'^2/2)- (24) 

Then the weights of the coordinate functions are W(yi) = 1 and VV(y2) = 4, while the weights 
of the coordinate vector fields are VV(g^) = 1 and = 4. 

Finally, let us write the Taylor expansion of the two vector fields fi,fo. 

d 

dyi 

d 

= (y?-22/12/2+ 2/i + o(|(yi, 2 / 2 )P) ^ 

We can see that /i has already weight 1, while the only part of weight 2 in /q is 2/i We can 
therefore conclude that the nilpotent approximation of /o,/i is given by 


fi 

fo 


fo=yi-^ and fi 

dy2 


d 

dyi' 


4.2 Order of the dilations 

We will analyze here the order of the dilations, that is the order of homogeneity of the volume 
form under the action of the dilations. This number will be crucial to find the order of degeneracy 
of the fundamental solution of the operator ©• 

Let us consider the dilations 5^. They were defined by introducing the notation x = 
(x^,..., x”^), where each component x* is a vector of length ki = di — dj_i. Then we set 
5e(x^, x^,..., x™') = (ex^, e^x^,..., e^x"*). Let N be the order of homogeneity of the volume 
form dxi A dx 2 A ... A dxn around the point xq, that is a number such that 

{6e)*{dxi A dx2 A ... A Xn) = dxi A dx2 A ... A dXn- 


Then N is given by 

m m 

N := ■ ki = (dimLj(xo) — dimLj_i(xo)). (25) 

i=l i=l 

Since this number is very important we give here some examples to understand its meaning. 

Example 2 (Continuation of Example [1]). As a first example we consider the one given in 
Example [T] We have already computed the filtration in the equations in (|23D , so we already 
know that the integers ki := dimLj — dimLj_i are 

ki = 1 k 2 = ks = 0 and k 4 , = 1. 

Therefore the order of the dilations isA^ = 1-1 + 4-1 = 5, as one can compute directly from 
the explicit expression of the dilations in (|24p . Notice that this number is much bigger then the 
dimension of the manifold. 

Example 3 (Sub-Riemannian manifold). Let us assume that the operator in ([T]) is induced 
by an equiregular sub-Riemannian manifold. In other words, we consider an operator without 
drift field, and the vector fields /i,...,/^ generate a completely non-holonomic equiregular 
distribution, A, of step m. Recall that the growth vector of the distribution is defined as the 
vector at any point q of the manifold given by 

(A(g), A‘^{q ),..., A^{q)) where A*+^ := [A, A*]. 
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Then the integers ki related to the filtration are the same defined by the growth vector, i.e. 
ki = dim(Aj) — dim(Aj_i), and the number N is exactly the homogeneous dimension Q of the 
manifold. More explicitly 

m 

N = Q = \ ■ k + 2 ■ k 2 + ■■■ rn ■ km = i(dim(Aj) — dim(Aj_i)). 

i=l 


In particular for the n+1 dimensional Heisenberg group the homogeneous dimension is Q = n+2. 

Example 4 (Linear case). As a last example, we consider an involutive distribution P, spanned 
locally by k constant vector fields, and assume that the drift field is linear in a neighborhood of 
xq. Without loss of generality we can assume that 


h 


d 

—— VI < i < A: and /o 

UXi 




d 

'dxj 


for some constants Aij. Under these assumptions, the only Lie brackets different from zero are 
the one involving only one vector field of the distribution and the drift field. Let us call A the 
n X n matrix with entries equal to Aij and B the nxk matrix that is the identity in the first k 
rows and is equal to zero in the last n — k rows. Then Hormander’s condition of hypoellipticity 
([^ becomes Kalman’s condition of controllability for linear control systems, that is the following 
condition on the rank of Kalman’s n x (nk) matrix 

T:ank[B,AB,A^B,...,A'^-'^B] = n. (26) 


The filtration is then completely determined and we have 


L 2 i-iixo) = L 2 i{xo) = span{A^B : 0 < j < i - 1}. 

Consequently the numbers d 2 i-i = d 2 i are determined by the rank of the Kalman’s matrix in 
(j26l) . where we stop the series of matrices at A^~^B. The numbers kj are zero if j is even, while 
if j = 2i — 1 they are the number of new linearly independent columns obtained by adding the 
matrix A^~^B to the previous one. The step of the distribution is then an odd number 2m — 1 
and N is equal to an odd sum of integers: 

m 

N = ^(2i - l)A: 2 j-i = 1 • A:i + 3 • ^3 + 5 • ^5 H-h (27fi - l)/c 2 m-i- 

i=l 

This number is already appeared in literature, namely in Chapter 4 of [T], where Agrachev, 
Barilari and Rizzi compute the small time asymptotics of the cost functional associated to any 
ample, equiregular geodesic 7 . In particular, if one computes the asymptotic of the cost func¬ 
tional associated to the ample, equiregular geodesics 'y{t) = xq, then the trace of the principal 
term of the expansion is exactly the same term N. We refer to |T] for an exhaustive presentation. 

5 Small time asymptotic on the diagonal 

We come back now to the perturbative method explained in Section [2l Let 

k 

^ - /o(v’) - ^ E ^ X 
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be the differential operator ([T]) on M+ x M and assume that fo, fi, ■ ■ ■, fk satisfy Hormander 
condition Q. Fix the graded structure, {U,x,w), around xq, induced by the filtration C in¬ 
troduced in Section [3] and let p{t, x, y) be the fundamental solution of the partial differential 
equation corresponding to a volume form /r, such that y = dx\ A ... A dxn in U. By Remark 12.21 
this assumption on y is not restrictive for our study. For 0 < e < 1 we dilate the space around xq 
with the dilation defined in Definition 13.61 then by Proposition 12.41 the fundamental solution 
in U of the operator 

+ (27) 


1 . 

dt 


— e 


2 = 1 


IS 


q^{t,x,y) = e^p{e^t,S^x,d^y) 


Vx, y 


(28) 


where N is the order of homogeneity of the volume form p computed in 

The fields fo, fi, ■ ■ ■, fk can be written in terms of the nilpotent approximation as 


fo — fo + 9o, 


fi = fi + 9i l<i<k, 


where qq and gi are vector fields of order less then or equal to 1 and 0 respectively. Then the 
operator in (I27p can be decomposed in a principal part, perturbed by a small operator: 




2=1 


d 


1 . 1 


^ ;2/o + oi^) + ^ E (7^* + ""( 7 ^ 


i=l 




2=1 

= ^ Lo-Xe 
at 


where is an operator that goes to zero as e and Lq is the operator defined by the nilpotent 
approximation: 

Lq ■■= fo + ■ 

2=1 

To this kind of operator we can apply Duhamel’s formula (1191) . but to this end we have to 
prove that there exists the fundamental solution of the principal operator d/dt = Lq. As we will 
prove now, this will follow by the property of hypoellipticity of the original operator, that are 
preserved by the nilpotent approximation, that defines Lq. The same statement can be found 
also in the paper by Bianchini and Stefani [7]. 

Proposition 5.1. Let /o,/i,...,/fc be the nilpotent approximation of the fields fo,fi,---,fk 
defined in Definition\4.1\ Then 


(i) for every bracket A such that ||A|| = i, then (Aj—A^) E Lj_i(xo) and A^^xq) = 0 whenever 
Aj(xo) E Tj_i(xo), where ||A|| denotes the weight of the bracket A defined in Definition 

E3 


(a) Assume Lm{xo) = Lie^^J/o, /i,..., fk}, then 

Liea,o{/o) fi, ■ ■ ■, fk} — Lie^ol/o) fi, ■ ■ ■, fk}- 


(29) 
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Proof. Let us prove the first statement. Let (JJ, x) be coordinates around xq adapted to the 
hltration {Lijj. Then for every i, 


Li{xo) = span | 


d 


XQ 


’ dxdi 


(30) 


XQ 


Let A be a bracket such that ||A|| = i, then as proved in |7j Theorem 3.1, 0{Af) < li, where O 
is the graded order associated to the graded structure induced by the filtration. Therefore there 


exist constants aj such that 






A. 

dxi 


(31) 


XO 


Notice that if two vector fields hi, /12 are homogeneous of graded order respectively rii and n 2 , 
then their Lie bracket is either zero or homogeneous of order ni + 712 . Then the Lie bracket A^ 
is either zero or homogeneous of order ||A|| = i. Therefore by equation (1311) . we have 


A/(*»> = E ^ 

j-.ij=ii 


(32) 


XQ 


By subtracting (I32p to m, we find that (Aj — A^) £ Lj_i(xo), because Li(xo) is obtained as 
in (f30]l . 

For the second statement, let (U, x) be as before and let V{xo) £ Lie^ol/oj /ij • • • j fk}- Then 
V{xq) = A^(xo), for some bracket A with ||A|| = j equal to the graded order of ^(xo). Then 
by expression (l30)l . there exist a* such that t^(a:o) = Since Lj{xo) C Lm(xo) = 

Lie^oi/o, /i, • • •, fk}, we have that V{xo) £ Liea;o{/o, /i, • • •, fk}- 

We prove the other inclusion by proving that Lj(xo) = span{A^(xo) : ||A|| < i}, for every i. 
We prove it by induction on i. 

For i = 1, Li{xo) = span{Aj(3:o) : ||A|| < 1}. Let Aj(xo) £ Li(xo), then by statement (z), 
{Af — A^){xo) £ Lo(xo) = {0}. Then Af(xo) = A^(xo) and the statement is true for z = 1. 

Assume that the statement is true for z — 1, then Lj(xo) = span{Aj(xo) : ||A|| < z} and 
(Af — A^)(xo) £ Lj_i(xo). By the induction hypothesis, there exists g £ span{A^-(xo) : ||A|| < 
z — 1} such that 

A/(xo) = A^-(xo) +g. 

And the statement is proved also for z. 

We conclude, since Lie^^J/o, fi,..., fk} = Lm(xo) C Lie^^ifo, /i,..., fk}- □ 

Corollary 5.2. The operator d/dt = Lq is hypoelliptic on M"'. 

Proof. By Hormander’s condition of hypoellipticity we know that hiex^ifo, fi, - - -, fk} = 
Then the hypothesis of statement (zz) of the Proposition are fulfilled and then also the nilpotent 
approximation is Lie bracket generating. Moreover, since /o(xo) = 0 = foi^o), the approxima¬ 
tion of the drift field can give some contribution in the generating process at xq only trough its 
Lie brackets with some other vector field. Then even Hormander’s condition of hypoellipticity 
(j^ holds at Xq for the nilpotent approximation, that is 

- /o, A, ■ ■ ■, A} = 

Using the lower semi-continuity of the rank, we can then find a small neighborhood U of xq 
where Hormander condition holds at any point. 
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Now by the homogeneity of the approximating system we know that 


Se*fo = e^/o and = efi VI < i < k. 


Therefore, since the differential operator commutes with the Lie brackets, we can extend Hor- 
mander condition, which holds on a neighborhood of xq, to the whole Euclidean space M”’ and 
the operator d/dt = Lq is hypoelliptic on M”. □ 


Remark 5.3. The assumption of fo{xQ) = 0 is necessaryfor the proof of this corollary. Indeed 
if fo{xo) 7 ^ 0 it could be that L is hypoelliptic, but Lq is not. 

For example, on the fields 


fi 


d 

dxi 


fo = (1 + Xi) 


A 

dxo 


satisfy Hormander condition Q, but this fails for their nilpotent approximation 


fi 


d 

dxi 


fo = 


d 

dx2' 


We can conclude that the principal part of the operator in (1271) admits a well defined heat 
kernel, qQ{t,x,y), for small time, that is given by the density function of the solution, of 
the stochastic differential equation in Stratonovich form 


k 

= foiOdt + ^ /i(0 o dwi 

i=l 


C(0) = X 


(33) 


where Wi is a 1-dim Brownian motion for every 1 < i < k. We can then apply the procedure 
introduced in Section [2] and we conclude by giving the asymptotic on the diagonal in xq of the 
fundamental solution p{t, x, y). 

Theorem 5.4. Let fo, fi, ■ ■ ■, fk be smooth vector fields that satisfy Hormander condition and 
{Lf, X, w) he a graded coordinate neighborhood induced by the filtration C defined in (I22p around 
a point xo, where /o(xo) = 0. Let qo{t,x,y) be the probability density function of the solution 
to equation (IMD with initial condition ^( 0 ) = x and assume that qo{l,xo,Xo) is strictly positive. 
Then the short time asymptotic on the diagonal of the fundamental solution, p{t, x, y), of the 
heat operator 

|-/o-^EV (34) 

2=1 

is given by 

p{t, Xo, Xo) = (1 + o(l)), (35) 

where N is the degree of homogeneity of the volume form dxi A ... A dxn under the action of the 
dilations 5^ computed in (I25p . 

Proof. By Duhamel’s formula (1191) the asymptotic of the fundamental solution q^ to the dilated 
operator as e tends to 0 is 


9e(L X, y) = qo{t, X, y) + qe* d4go(T x, y), (36) 

provided that the remainder term q^, * X^qo is negligible for e small. 
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Let Z be the operator We have to prove that 



s,XQ,y)ZqQ{\ 


s,y,xo)dyds 


is bounded uniformly in e. We split the integral between the integration on 0 < s < ^ and the 
integration on ^ < s < 1. For the first integral, since we are integrating on a bounded set where 
1 — s is far from zero, ZqQ{l — s,y,XQ) is uniformly bounded by a constant C. Then 



qe{s,xo,y)Zqo{l 


s,y,xo)dyds < 


< 



qe{s,xo,y)dyds 

e^p{e^s,xo,5ey)dyds 


where the last inequality follows since p has integral equal to 1 over the whole space. 

For the other part of the integral, the smooth function (7e(s,Xo,y) is uniformly bounded on 
U, since s is far from 0. Moreover, the operator Z is a. combination of polynomials, that are 
bounded on the bounded set U and of derivations of go of ^-^y order. So let us control the 
integral 



s,y,xo)dyds 


for every 1 < i < n. By an integration by parts, it is equal to 


/ / qoi'i--s,yui,XQ)-qo{l-s,yu2,xo)dyi---dyi---dynds 

Ji/2 Jc/nR"-! 


where yjj^ and yu 2 a-re the two extremals of the set U on the line of integration in yj, while dyi 
denotes that we are not integrating in dy^. This integral is bounded since go has integral equal 
to 1 on the whole space. 

For the higher derivatives, let us assume by induction that for every set of indexes I, with 
l-f| = j) we have proved that the integral 



s,y,xo)dyds 


is bounded. Let J be a set of indexes such that | J| = j + 1, and that it differs from I by an 
index i, then we can integrate by parts in the coordinate xt the whole integral 


f f qe{s,xo,y)-^qo{l - s,y,xo)dyds = 

J 1/2 Ju 9 xj 

= [ f qe{s,xo,y) -^qo{l - s,y,Xo) 

Ji/2 JunRn-i 


VUi 


dyi---dyi--- dynds- 


yu2 


f f [ 

Ji /2 JunM.n -1 Ju 


_i JunR ^xi 


dq^ 

xq, y)—go(l - s, y, xo)dyds 


and this is bounded by the induction hypothesis and since g^ and ^ are bounded away from 

s = 0. 
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Now let us come back to the asymptotic in (I36p . By the definition of q^, if in (|28p we fix 
t = 1, X = y = xq, and we let e go to zero as e = \/t, we notice that 


p(t,xo,xo) = 


q^(l,xo,xo) 


Vt 


■N 


The desired small time asymptotic on the diagonal is then determined by the asymptotic 
for and we find 

p{t, Xo, Xo) = --(1 + o(l)), 

which is well defined since by hypothesis the leading term g'o(l)Xo,xo) doesn’t vanish. □ 

Remark 5.5. The conclusion of Theorem 15.41 holds also for the operator 

4 - /o - J '^ifi + (div;,/i)/i) 


(37) 


2=1 


where p is any volume form on M and fo, fi, ■ ■ ■, fk satisfy the assumptions of Theorem 15.41 
Indeed we can treat the operator in (1371) as an operator of the form 






2=1 


where we take /o = /o + ^ X)i)=i(div^/j)/j, but in this case we could have /o(xo) ^ 0. How¬ 
ever, if fo, fi, ■ ■ ■, fk satisfy Hormander condition ([2]), then it is satisfied also by fo, fi, ■ ■ ■, fk- 
Moreover, the filtration C generated by fo, fi, ■ ■ ■, fk is equal to the filtration C generated by 
fo, fi, ■ ■ ■, fk, since fo differs from fo by a linear combination of the vector fields fi,..., fk- Then 
the approximating system fo, fi, - - -, fk is the same as the approximating system fo, fi, - - -, fk, 
hence the operators in (I37p and in ()34p have the same principal part. The proof of Theorem 15.41 
for the operator in (|37p now follows as the one for the operator (I34D . 


6 The principal operator and the associated control system 

In this section we are going to investigate the conditions for the positivity of the heat kernel, 
qo{t,x,y), that we have introduced in the last section. 

Let fo, fi, - - -, fk satisfy Hormander condition ([2]) and consider the principal operator 

2=1 

defined by the approximating system of the original system of vector fields. As already pointed 
out, it admits a smooth fundamental solution given by the probability density, qo{t,x,y), of the 
process to be at time t in the point y starting from the point x, where is the solution of 
the stochastic differential equation 


k 

= fo{it)dt + ^ fi{it) o dwi- 

2 = 1 


(39) 
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In their famous work m Stroock and Varadhan characterized the support of qo{t,x,y) and 
they showed that it is the set of reachable points from x of the following associated control 
problem: 

k 

X = foix) + ^ Uifi{x) ( 40 ) 

i=l 

where x : [0, t] M”' is a curve in M"' and u = {ui ,..., u^) € t];are bounded controls. 

Unfortunately, Stroock and Varadhan’s result holds only for globally bounded vector fields, 
with bounded derivatives of any order. Since our vector fields are polynomial, they don’t satisfy 
such assumptions and we can not directly apply the result of the support theorem. We will see 
in a moment how we can adapt their procedure to our system, but hrst we introduce two simple 
lemmas that will simplify our study. 

Lemma 6.1. Let qo(t,x,y) be the probability density of the solution to equation (I39p . Then for 
every e > 0, for every t > 0 and for all x € ML it holds 

qo{t, xq, x) = e^qo (e^t, xq, d^x) . 

In particular, we have the following equivalence 

qo{l,Xo,Xo) > 0 qo{t,Xo,Xo) Vt > 0 . ( 41 ) 

Proof. This is a corollary of Proposition 12.4[ since by definition of /oj /ij • • •, /fc we have 

= Lq. 

€ 

□ 


Lemma 6.2. Consider the control problem ()40]) on M”. Let yi,y 2 G K"' and T > 0 be fixed 
and assume there exists a curve y : [0,T] —)• M” that satisfies the control problem (I40D for some 
control function u G and such that y(0) = yi and y{T) = y 2 . 

Then for any M > 0 the curve x{t) := 6 m {u (]\p')) *-5 an admissible curve for the control 
problem defined on [0, M^T] with control u{t) := -^u (t/M'^'j, that connects xi := dMivi) with 
X 2 ■= 6 M{y 2 ) in time M'^T. 

In particular, if the control problem ini) is controllable in a neighborhood U of xq in time 
T, then it is controllable in 6 m{U) in time M‘^T. 


Proof. The boundary conditions are easily satisfied, since x{0) = 6Miy{0)) = 6 M{yi) = xi and 
x{M‘^T) = 6 M{y{T)) = 6M{y2) = X2. Moreover, x{t) is an admissible curve for the control 
problem with control -^Ui {t/M'^^, indeed by the homogeneity of the approximating system we 
have 


x{t) = 


1 

M2 

1 

M2 

1 

M2 


6m* [y {t/M"^)) 


6m* 


fo {y {t/M"^)) {t/M"^) fi {y {t/M"^)) 


i=l 

k 


M2/o {y it/M^)) +Y,Mui {t/M^) h {y it/M^)) 


2=1 


/o {x{t)) + ^ —Ui {t/M"^) fi {x{t)). 


2=1 


□ 
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We generalize here Stroock and Varadhan’s support theorem for the control problem defined 
by the approximating system. Indeed their theorem in |13) can be applied only to fields that 
are bounded and Lipschitz, with bounded derivatives of the first and second order, and does 
not hold for general unbounded fields. In our case, the helds are polynomial, so in general not 
globally bounded. However they are still very particular, because, since /o has weight 2 and 

fi has weight 1, for every 1 < i < k, every component of these fields depends only on the 

"(i) 

coordinates of less weight, i.e. f> does not depend on ... ,Xn- This will allow us to 

modify the proof of the support theorem in order to extend it to our case. 

Definition 6.3. Consider the control problem (14011 and let us call Xu{t) the solution correspond¬ 
ing to a control u. The reachable set of the control problem in time t from x is the set, 

At{x) := I?/ G M” : 3u G L°°([0, t];M^) such that Xu(0) = x and Xu{t) = ?/| . 

Lemma 6.4. Let Xq,Xi, ..., be smooth vector fields on M"', that satisfy Hormander condi- 
tion, and such that every j-th component Xf of Xi, for 0 < i < k, does not depend on the 
coordinates Xj,... ,Xn, but only on the first coordinates xi,... ,Xj-i. Let be the solution of 
the stochastic differential equation 


k 

dft = M^t)dt + J2x^i^t)dwt 

i=l 


where Xq stands for the vector field whose j-component is given by 


-U) 


X-> = X, 


(i) 


+ 


i+l 1=1 


(0 


dX, 


U) 


dxi 


Let p{t,x,y) be the probability density of to be in y at time t starting from the point x. Let 
At{x) be the reachable set at time t from x of the associated control problem 

k 

X = Xo{x) + ^ Ui{t)Xi{x) (42) 

i=l 

where u = (ui, ... ,U}f) G L°°([0, t];M^) is a control function. Then 

supp{p{t,x, •)) = At{x). 


Proof. Stroock and Varadhan have proved this theorem under the assumption that the fields 
are Lipschitz and globally bounded, together with their derivatives of first and second order. 
Following their proof in |13) . we have to show that for a dense set of controls u and Ve > 0 


PxiW^t - xt\\ < e) := P(||6 - xt\\ < e| ^0 = a;) > 0, 


where xt is the solution of (|42p starting at x. In particular let us take G M^), with 

"0(0) = 0, and let Xt be the solution of ()^2l) starting at x with control uft) := fiift). Then for 
all e > 0 we show that 


PxiWft - xtW < e I \\wt- iftW < 5) 1 


(43) 


as 5 0. Indeed 


Px (116 - xt\\ < e) = (116 - xt\\ < e I \\wt - 6*11 < <5) ■ P{\\wt - 6*11 < ^) 
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and P (ll^t — V’tll < <5) > 0 for every 5 > 0. 

Strooch and Varadhan proved (j33|) under the boundedness assumption that we don’t have 
directly, but we will recover it by iterating a conditional probability. Indeed, notice that by our 
assumption the first component of every vector field, , does not depend on any coordinate, 
so they are actually constant and they trivially satisfy Stroock and Varadhan’s assumptions. 
Then the limit in (|43p holds for the process ~ 

Moreover let us assume, by induction, that the first j — 1 components of live in a bounded 
set. Then the components Xf are Lipschitz and bounded, together with their derivatives of 
any order, and we can apply Stroock and Varadhan’s theorem to the j-th component of then 

Px < e I \\wt - -iptW < 5, < e VI < / < j) ^ 1 


as 5 \ 0 . 

The proof of (j43|) now follows using an iterated conditional probability, indeed in general for 
every measurable set Mi,..., An, B, it holds 





Then 


\i--i 

\j.3 


j=i \ 


.BnMij P{Ai\B) 

.BnMinMa I P {A 2 \b r\ Ai) P {Ai\B) 


i -1 ' 

Bnf]Ai 

1=1 


P{\\Ct - xtW < e I \\wy - iptW < 6) = 

n 

= < e I Wwt-iitW < 5, -x^\\ < e VI < Z < j) ^ 1 

i=i 


as 5 0 and we have proved f)^3ll in our case. □ 

We are now ready to show a condition for the positivity of the fundamental solution of the 
approximating differential operator. 

Theorem 6.5. Let qQ{t,x,y) be the fundamental solution of (1381) . If the reaehable set At{xo) of 
the assoeiated eontrol problem ini) is a neighborhood of xq for some t > 0, then 50(15 xq, xq) > 0. 

Proof. By Lemma l6.ll it is enough to prove that qo(T, xq, xq) > 0 for some T > 0. We will 
choose T = 2t. Moreover, by Lemma 16.21 if Vli(xo) is a neighborhood of xq for some f > 0, it is 
a neighborhood for every t > 0. 

Assume by contradiction that qQ{2t,XQ,XQ) = 0. By Chapman-Kolmogorov equation we 
know that 

0 = qo{2t,xo,xo) = / qo{t,xo,y)qoit,y,xo)dy = /_ qo{t,xo,y)qo{t,y,xo)dy, 

-Ik" J At(xo) 
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where we can restrict the space of integration, since by Lemma 16.41 supp{qo{t,Xo, •)) = At{xo). 
Then for all y S At{xo) we have qQ{t,y,Xo) = 0. The function q(t,x,y) := qo(t,y,x) is the 
fundamental solution of the adjoint operator to Lq, that is 



2 = 1 


Then q{t, x, y) is the probability density function of the stochastic process solution of the 
stochastic equation 


k 

dlt = -fo{it)dt + ^ fiiit) o dwi{t) 

2=1 

io = X. 

By contradiction we have assumed that q{2t,XQ,XQ) = 0, then again by Chapman-Kolmogorov 
equation we have 

0 = q{2t,Xo,xo) = I q{t,Xo,y)q{t,y,xo)dy 
J At (xo) 

where At{xQ) is the reachable set in time t from the point xq of the associated control problem 

k 

X = - fo{x) + '^Ui{t)fi{x). 

2 = 1 

It follows that qo{t,Xo,y) = q{t,y,Xo) = 0 for all y £ At{xo). Since the point xq is a stationary 
point for the control problem, then xq £ ^t(xo), for every t > 0. By Krener theorem, xq is in 
the closure of int(^ 4 (xo)), then At{xo) H At{xo) has non zero measure and qo{t,Xo,z) = 0 for 
all z in this intersection. This is a contradiction to the support theorem. 

We conclude that, if is a neighborhood of xq, then qo(t, xq, xq) > 0 for all t > 0. □ 

Remark 6.6. In view of Theorem 15.41 and Theorem 16.51 we can conclude the following properties 
about the asymptotic of the fundamental solution p of the operator ©• 

(i) if the control problem associated with the original system: 

k 

X = foix) + '^ Ui{t)fi{x) (44) 

2=1 

is not controllable around xq, that is At{xQ) is not a neighborhood of xq, then 

p{t,xo,xo) = 0 Vt > 0. 

Indeed, even if the support theorem can not be applied to this system, we still have that 
supp(p(t, Xo, •)) " 44 ( 2 ^ 0 ), therefore xq is on the boundary or out of the support. Since 

p{t,XQ, ■) is smooth, the conclusion follows. 

(ii) if the control problem (j44p is controllable in xq, then we study the controllability of its 
nilpotent approximation, defined by the fields /o, /i, • • •, /fc introduced in Sectional 
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(11.1) if the approximating control problem 

k 

X = fo{x) + ^ Ui{t)fi{x) (45) 

is controllable around xq, then the asymptotic is given in Theorem 15.4[ where we 
see that the fundamental solution on the diagonal in xq blows up for small t as the 
rational polynomial for a constant oq depending on the chosen volume and 

on the approximating system (|45p . The order N is determined by the Lie algebra 
generated by the fields /o, /i, • • •, /fc at xq as explained in formula (I25p . 

(11.2) it the approximating control problem (I45p is not controllable, then either p{t, xq,Xo) = 
0 for every t > 0 (but this is not possible if the fields fo, fi, ■ ■ ■, fk are bounded, with 
bounded derivatives, as explained in the proof of Theorem 16.511 . or p{t, xq,xq) goes to 
infinity for small t faster than as shown in the following Proposition. In [5] the 
authors show an example, where the asymptotic goes to infinity even exponentially 
fast. 

Proposition 6.7. Assume that the eontrol problem (j43|) is controllable around xq, but the ap¬ 
proximating control problem (05]) is not. If p{t,xo,xo) > 0, then limt\ot^p{t,xo,xo) = +oo, 
where N is defined in (1251). 

Proof. To show this we need to prove that all the coefficients in the asymptotic (13511 vanish. The 
asymptotic was found introducing the fundamental solutions q^: and qq and by using Duhamel’s 
formula (|19p . By iterating it we can achieve a better approximation of the asymptotic and find 
the higher coefficients. Indeed,we obtain 

3 

Qt = qo + '^qo {*^eqoy + qe i*A:^qoy^^ 

i=l 

for every j S N, where means that we iterate the convolution i times. We have to show 

that all the terms go {*^eQoy vanish at the point (l,3:o,xo) for every z E N. 

By Lemma 16.41 we know that qQ{t,xo,y) = 0 for every y E At{xo) and all t > 0. Let us 
assume by induction that 

qo{*^eqoy {t,xo,y) = 0 '^y e Afixof, '^t > 0. (46) 

- C 

Then for every y E At{xo) we have 

go (*T’ego)*+^ {t, xq, y)= /_ go (*T’,go)* (s, xq, z) X^qo{t - s, z, y)dzds, 

Jo JAs{xo) 

where the integral can be computed just on As{xo) by the induction hypothesis. But go(t — 
s,z,y) = 0 on .As(xo) by Chapman-Kolmogorov equation, then also all its derivatives vanish 
there. Then the integral is zero and we have proved (I46p for every i. 

Since go is smooth and xq is on the boundary of Afixo) , then equation (14611 holds also in 
the point {t,Xo,Xo) for every t > 0, that means that all the coefficients of the asymptotic (I35p 
are zero. □ 
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7 Examples 


We end this paper with a study of some known examples, to understand better the meaning of 
Theorem 15.41 

Example 5 (Continuation of Example [T]). We complete the study of Example [T] We have 
already computed the principal part of the dilated operator L^, that is 

dt ^^dx2^2dxl' 

This operator is indeed hypoelliptic since /i = and [/i,[/i,/o]] = span the whole 

tangent space in every point. Let qQ(t,x,y) be the density function of the solution ^{t) of the 
stochastic equation f)33p . that in coordinates is given by 

= dwi, d ^2 = x\dt. 


We can see that the second coordinate is actually deterministic and has positive derivative, so 
it can only increase. Consequently, if a path starts in xq = (0,0) the solution ^(t) will almost 
surely never come back to xq again, indeed xq is on the boundary of the support of qo{t,Xo, ■). 
Then the hypothesis of Theorem 15.41 that requires go(l) 3^0) 3^o) > 0 is not fullfilled. 

Since the original control problem is controllable, this is an example of the type (ii.2). 


Example 6 (Sub-Riemannian manifold: continuation of Example [3|). The study of the asymp¬ 
totic on the diagonal of the heat kernel on 3D contact sub-Riemannian manifolds has already 
been performed by Barilari in [2]. The nilpotent approximation of a sub-Riemannian manifold 
of dimension 3 is isometric to the Heisenberg group. Let us represent the Heisenberg group as 
with coordinates {x,y,z), then the approximating system can be written as 


d y d 


d X d 
dy 2 dz 


As one can easily verify, the order of homogeneity of the volume form is given by = 4, as 
computed also with the general formula (f25l) . 

The principal part of the dilated operator is hypoelliptic and symmetric, the associated 
control problem is then controllable, so there exists a well-defined symmetric heat kernel, that 
is positive for every t > 0 as seen in Theorem 16.51 The hypothesis of Theorem 15.41 are then 
fullfilled and we find that the asymptotic on the diagonal of the original heat kernel p{t, q, q') 
has the following order: 


p{t,q,q) 


ao{q) +o(l) 


for a smooth function ao(<?) > 0 on the manifold. 


This was the same order found by Barilari in |2]. 
This example is of the type (ii.l). 


Example 7 (Ben Arous and Leandre). We consider here an example studied by Ben Arous and 
Leandre in |5]. Consider the space with coordinates {x\,X 2 ) and let 


k=x 


-A 

^ dx2 ’ 


h 


d 

dxi ’ 


f2 = X 


^dx2 ’ 


where a and h are positive integers. Then xq = (0,X2) is a stationary point of the drift field for 
any X 2 S K. The operator L = Jq + ^(/i + /|) satisfies even the strong Hormander condition, 
i.e. the fields /i ,/2 alone are Lie bracket generating. 

In [6], it is shown the complete behaviour of the heat kernel p(t, x, y) of this operator on the 
diagonal. We summarize here the most interesting properties for our study: 
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Theorem 7.1. 1. If h < a + 1, then there exists a constant K(a, b) > 0 such that 


p{t,xo,xo) ~ 


K{a, b) 


( 47 ) 


2. If b> a+1 and a is even, then the fundamental solution p(t, xq, xq) decreases with expo¬ 
nential velocity. 

The results found in this paper agree with the statement of the theorem. Indeed, let us 
consider the filtration C given by the vector fields /o,/i,...,/fc at a point xq = (0,X 2 ). Let 
m := min{a + 2 ; 6 + 1 }. It is easy to verify that the subspaces of the filtration C are 


Li{xo) 


M X {0} if 1 < i < m 
if f = m. 


Accordingly the coordinate xi has weight 1, while the coordinate X 2 has weight m and the order 
of homogeneity of the volume form is given by 


f 6 + 2 if 6 <a + 
\a + 3 if 6 >a + l 


To determine the nilpotent approximation, it is convenient to divide the study in 3 cases, de¬ 
pending on the value of a w.r.t. 6. 

If 6 < a + 1, then m = 6 + 1. The nilpotent approximation is obtained by taking the Taylor 
expansion of the field /o of order 2 and the Taylor expansion of order 1 of /i and / 2 . Then we find 
that /i = /i, /2 = /2 and /q = 0. The principal part of the dilated operator is ^(/i + f^), that 
is hypoelliptic. Then there exists a well-dehned heat kernel, qQ{t,x,y), and, since the associated 
control system is controllable, qQ{t,XQ,Xo) > 0 for every t > 0. Then the hypothesis of Theorem 
lOl are fulhlled and we find that the small time asymptotic of the fundamental solution of L has 
order N/2 = that is exactly the one given in (1471) . 

If 6 = a + 1, then m = 6 + 1 and the nilpotent approximation is equal to the fields fi, f 2 , fo 
themselves. Since /i ,/2 are Lie bracket generating the associated control system 


x = fo + Ulfl + Uxf2 

is still controllable, then the heat kernel qo{t,Xo,Xo) is positive for every t > 0 , and we obtain 
again the statement (SZl). 

If 6 > 0 + 1 , then /i = /i, /o = fo and /2 = 0. The principal operator /o + 5/1 is still 
hypoelliptic, but if a is even, then the heat kernel qo is zero in xq = (0, X 2 ) for any t > 0. This 
is because a.e. path starting from xq will never come back to xq again, since the drift /q makes 
the first coordinate increase, if xi becomes different from 0. Then we can not apply Theorem 
EM and indeed Ben Arons and Leandre have shown an exponential decrease in this case. 
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